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Abstract— The use of Space-Time Block Codes (STBCs) in- 
creases significantly the optimal detection complexity at the 
receiver unless the low-complexity decodability property is taken 
into consideration in the STBC design. In this paper we propose a 
new low-complexity decodable rate-5/4 fuU-diversity 4x4 STBC. 
We provide an analytical proof that the proposed code has 
the Non-Vanishing-Determinant (NVD) property, a property that 
can be exploited through the use of adaptive modulation which 
changes the transmission rate according to the wireless channel 
quality. We compare the proposed code to the best existing low- 
complexity decodable rate-5/4 full-diversity 4x4 STBC in terms 
of performance over quasi-static Rayleigh fading channels, worst- 
case complexity, average complexity, and Peak-to-Average Power 
Ratio (PAPR). Our code is found to provide better performance, 
lower average decoding complexity, and lower PAPR at the 
expense of a slight increase in worst-case decoding complexity. 

Index Terms — Space-time block codes, low-complexity decod- 
able codes, conditional detection, nonvanishing determinants. 

I. Introduction 

Space-time coding techniques have become common-place 
in wireless communications standards [1], [2] as they provide 
an effective way to mitigate the fading phenomena inherent 
in wireless channels. However, the use of Space-Time Block 
Codes (STBCs) increases significantly the optimal decoding 
complexity at the receiver unless the low-complexity de- 
codability property is taken into consideration in the STBC 
design. We distinguish between two decoding complexity 
measures, the worst-case decoding complexity measure and 
the average decoding complexity measure. The worst-case 
decoding complexity is defined as the minimum number 
of times an exhaustive search decoder has to compute the 
Maximum Likelihood (ML) metric to optimally estimate the 
transmitted symbols codeword f3], fTl, or equivalently the 
number of leaf nodes of the search tree if a sphere decoder is 
employed, whereas the average decoding complexity measure 
may be numerically evaluated as the average number of visited 
nodes by a sphere decoder in order to optimally estimate 
the transmitted symbols codeword JS]. It is well known that 
Complex Orthogonal Design (COD) codes f6l, |T| are Single- 
Symbol Decodable (SSD) for general constellations. But if one 
considers the case of rectangular QAM constellations, the de- 
tection of each complex symbol reduces to separate detection 
of two real symbols which can be effectively performed via 
two threshold detectors (or equivalently PAM slicers) without 



any decoding complexity. On the other hand, the rate of 
square COD codes decreases exponentially with the number 
of transmit antennas I?), which makes them more suitable for 
low-rate communications. In this paper, we address the issue of 
increasing the rate of the 4x4 COD code at the expense of an 
increase in decoding complexity while preserving the coding 
gain and the Non- Vanishing-Determinant (NVD) property |8|. 
This property can be exploited through the use of adaptive 
modulation which varies the transmission rate (through the 
choice of the modulation order) according to the wireless 
channel quality. We propose a new low-complexity decodable 
rate-5/4 full-diversity 4x4 code that encloses the rate-3/4 
COD code in f7\ and we provide an analytical proof that the 
proposed code has the NVD property. 

We consider the case of square QAM constellations and 
compare the proposed code to the best existing low-complexity 
decodable rate-5/4 full-diversity 4x4 STBCs in terms of 
performance over quasi-static Rayleigh fading channels, worst- 
case decoding complexity, average decoding complexity, and 
Peak-to-Average Power Ratio (PAPR). Our code is found to 
provide better performance, lower average decoding complex- 
ity and lower PAPR at the expense of a slight increase in 
worst-case decoding complexity that will only penalize the 
proposed code for high-order constellations. 

The paper is organized as follows: In Section II, preliminar- 
ies on COD are provided. In Section III we introduce the new 
code and derive its low-complexity decodability property. In 
Section IV, performance comparisons by means of computer 
simulations are provided. We give our conclusions in Section 
V, and the proof of the NVD property for the proposed code 
is provided in the Appendix. 

Notations: 

Hereafter, small letters, bold small letters and bold capital 
letters will designate scalars, vectors and matrices, respec- 
tively. If A is a matrix, then A^, A-^ and t/ef[A] denote 
the hermitian, the transpose and the determinant of A, re- 
spectively. We define the vec{.) as the operator which, when 
applied to a m x n matrix, transforms it into a mn x 1 vector by 
simply concatenating vertically the columns of the correspond- 
ing matrix. The (g) operator is the Kronecker product and Skj is 
the kronecker delta. The sign{.) operator returns 1 if its scalar 
input is > and -1 otherwise. The round{ ) operator rounds its 



argument to the nearest integer. The and operators 
denote the real and imaginary parts, respectively, of their 
argument. The (.) operator when applied to a complex vector 
a returns [3? (a^) ,3(a^)] and when applied to complex 

matrix A returns [3? (A-^) , 3 (A-'")]"^. For any two integers 
a and b, a = b (mod n) means that a — & is a multiple of n. 

II. Preliminaries 
We define the MIMO channel input-output relation as: 

Y = X H + W (1) 

TxNr TxNtNtxNr TxNr 

where T is the number of channel uses, Nr is the number 
of receive antennas, Nt is the number of transmit antennas, 
Y is the received signal matrix, X is the code matrix, H is 
the channel matrix with entries hki ^ CAf{0, 1), and W is 
the noise matrix with entries Wij ~ CAf{Q, Nq). In the case 
of Linear Dispersion (LD) codes [9], an STBC that encodes 
2K real symbols is expressed as a linear combination of the 
transmitted symbols as: 



2K 



k=l 



(2) 



with Xk G M and the fc = 1, 2K ecre T x Nt complex 
matrices called dispersion or weight matrices that are required 
to be linearly independent over M. 

The MIMO channel model can then be expressed in a useful 
manner by using ^ as: 



Y 



2K 

E 

k = l 



(3) 



Applying the vec{.) operator to the above equation we obtain: 



2K 

vec{Y) = J2 (Ia'. ® 0k) vec (H) Xk + vec(W). (4) 

fc=i 

where Ijv^ is the x iV^ identity matrix. 
If Yi, h.i and designate the i'th column of the received 
signal matrix Y, the channel matrix H and the noise matrix 
W respectively, then equation Q can be written in matrix 
form as : 
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Thus we have: 



y = ns 



(5) 
(6) 



A real system of equations can be obtained by separating the 
real and imaginary parts of the above to obtain: 

(7) 



y 

2Nr-Txl 



n S- 

2Nr-Tx2K 



w 

2Nr-Txl 



n = [Qi Q2] 



(8) 



Assuming that N^T > K, the QR decomposition of "fi yields: 

TR] 


where Qi e R^^^^x^^,Q2 e r2N^tx(2N^t-2K) ^ qTq^ ^ 
I, i = 1, 2, R is a 2K x 2K real upper triangular matrix and 
is a {2NrT - 2K) x 2K null mati'ix. Accordingly, the ML 
estimate may be expressed as: 



— arg min||y — QiRs| 



(9) 



where C is the vector space spanned by information vector s. 
Noting that multiplying a column vector by a unitary matrix 
does not alter its norm, the above reduces to: 

Rs||2 



s^'' = arg min||y' 



(10) 



where y' = y. 

Complex Orthogonal Design (COD) codes: 

For COD codes, the weight matrices satisfy Q: 

Af^j +^f^, -2%I^, (11) 
In this case, if hi is the i'th column of the equivalent channel 
matrix H in (|6]l, then it is straightforward from ( fTT) to prove 
that: 

(12) 



h]^hi + hfhk = 0, k ^ I 



or equivalently: 



3? 



{fiffi/} = V A: 7^ /. 



(13) 



in terms of the columns of % namely fi, the above reduces to 

hlhi = oy k^i. (14) 

The orthogonality of the columns of 71 is inherited by the 
upper triangular matrix R in ( fTO] ) which becomes simply 
a diagonal matrix. In practical communication systems, the 
transmitted symbols are drawn from complex constellations 
and thus the code matrix X can be seen to encode K complex 
symbols Si where X2i-i and X2i are the corresponding real 
and imaginary parts respectively with i = 1,...,K. The 
choice of the used constellation plays a key role in deter- 
mining the worst-case decoding complexity, because if the 
complex symbols Sj are drawn from a general constellation, 
the corresponding real and imaginary parts cannot be detected 
independently and thus the worst-case decoding complexity 
is 0{M) where AI is the size of the constellation. On the 
other hand if the complex symbols Si are drawn from a 
rectangular QAM constellations, the ML decoding process 
of each complex symbol Si reduces to separate detection 
of the real and imaginary parts and the COD code can be 
decoded via 2K PAM slicers as shown in ( fTSj ). It is worth 
noting that the PAM slicer equations ([TSj require only a 
fixed number of simple arithmetic operations, which does not 
grow with the size of the rectangular QAM constellation, 
and therefore according to the definition of the worst-case 
decoding complexity they are considered of complexity 0(1). 



xf^ = sign {y'Jr.^i) x min |2 round{ {y'Jr,,, - 1) /2) + l|, VM - 1* 



Vi 1, 



.2K. 



(15) 



We assumed here that the underlying constellation is square such as 4-/16-QAM constellations, which is usually the case in practical scenarios 



^ne'vv(s) 



xi +jx2 - jxwe^"^ 
-X5 + jxa +jxge'^ 



X3 + jX4 

xi - 3x2 - jxloe 
e^-^ixr+jxs) 
X5 - jxe + jxge 



34> 



X5 +jX6 + jxge^''' 

e^^{-X7+jxs) 
x\ - jx2 +jxioe^'^ 

-~X3 + jX4 



-e^^{x7+jxs) 
-X5 - jxe +jxge^'' 

X3 + jX4 



(17) 



xi + jx2 + jxioe- 



III. The Proposed Code 

From the previous section, a judicious structure for a new 

high-rate, low-complexity code would be to embed a COD 

code into a new higher rate STBC. The resulting code will 

enjoy the low-complexity decodability through the use of 

conditional detection ifTOll . ifTTl . The proposed code denoted 

X.„g„ is expressed as: 

X„„.(s) =0(a;i, . . . ,X6)+ 

■A. (16) 
e^''' {112x7 + RaXs + Rixg + Rsxio) R4 

with s = [si,S2], si = [xi,...,xe], S2 = [x7,...,xiq] 

and (p is chosen to maximize the coding gain. The 

0{xi, . . . ,X2{a+i)) matrix denotes the square COD code 

for the case of 2" transmit antennas and is the matrix 

representations of the Clifford algebra generator 7^ with 1 < 

i < 2a + 1 Q. The proposed code matrix takes the form 

of ( [T7] l. It was verified through exhaustive search that taking 

(f> = ^cos~^(l/5) maximizes the coding gain and that it 

remains constant up to 64-QAM unnormalized constellations. 

An analytical proof that the coding gain of the proposed code 

is constant over unnormalized QAM constellations is provided 

in the Appendix. The corresponding upper-triangular matrix R 

takes the form below: 
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(18) 



where x indicates a possible non-zero position. The de- 
coder exploits the structure of the upper triangular matrix R 
by computing the ML estimates of the orthogonal symbols 
(xi, . . . , xq) assuming that a given value of (2:7, . . . , xiq) is 
transmitted. In the case of the square QAM constellations, 
the ML estimates of the orthogonal symbols {xi,...,xe) 
assuming the knowledge of {x7,...,xiq) can be obtained 
exactly as in ([T5| with the only difference that y'^ will be 
replaced by Zi = y[ — X]^=7 ^i.j^j- Therefore, the worst-case 
decoding complexity of the proposed code is equal to AP as 
the hard PAM slicers prune the last 6 levels of the original 10 
levels real valued search tree. 

IV. Numerical and Simulations Results 

Table]!] at the top of the next page, summarizes the compar- 
ison between the proposed code and the rate-5/4 punctured 



version of the code in |12| in terms of square QAM constel- 
lations worst-case complexity, the minimum determinanj^^nd 
the PAPR. The minimum determinant is defined as: 

Min det= min |t/er [(X(As))] I = (19) 

AseAC\{0} 

where As = s — s', AC is the vector space spanned by As, 5 
is the coding gain and the PAPR is defined as: 



PAPR„ 



max|X(t,n)|^ 



(20) 



r-iEE{|X(i,n)|2} 
t 

where i G {1,...T} and n e {l,...Nt}. Due to the 
symmetry between transmit antennas, the subscript n will 
be omitted. The proposed code has a higher coding gain, 
lower PAPR at the expense of a slight increase in worst-case 
detection complexity that will penalize our code for the 64- 
QAM constellation and above. 

Simulations are carried out in a quasi-static Rayleigh fading 
channel in the presence of AWGN for 2 receive antennas and 
4-/16-QAM constellations. The ML detection is performed 
via a depth-first tree traversal with infinite initial radius 
sphere decoder. The radius is updated whenever a leaf node 
is reached and sibling nodes are visited according to the 
simplified Schnorr-Euchner enumeration flSl. Fig.]T]illustrates 
the performance comparison in terms of Codeword Error 
Rate (CER) of the proposed code and the rate-5/4 punctured 
version of the code in fV2\. One can notice that the proposed 
code offers better performance especially at high SNR. This 
can be interpreted by the superiority of the coding gain of our 
code. Furthermore, as can be seen form Fig. ]2] our code can 
be decoded with lower average complexity. 

V. Conclusion 

In this paper, we have introduced a new low-complexity 
decodable rate-5/4 full-diversity 4x4 STBC that encloses 
the rate-3/4 COD and retains its coding gain. The coding 
gain is proved analytically to be constant over unnormalized 
QAM constellations. The proposed code is compared to the 
rate-5/4 punctured version of the low-complexity decodable 
4x4 STBC in O!). The proposed code showed better CER 
performance, lower PAPR, lower average complexity at the 
expense of a slight increase in worst-case detection complexity 
that affects our code for high-order QAM constellations. 

Appendix 

In the following, we will prove that choosing = 
i cos" ^(1/5) indeed maximizes the coding gain 5 (which is 
equal to IQ^) for unnormalized QAM constellations and thus 
guarantees the NVD property for the proposed code. 

*For consistency, the minimun determinant is evaluated at fixed average 
transmitted power per channel use for both codes 



TABLE I 

SUMMARY OF COMPARISON IN TERMS OF WORST-CASE COMPLEXITY, MiN det AND PAPR 



Code 


Worst-case 
Complexity 


Min det (= VS) 
for QAM constellations 


PAPR (dB) 


QPSK 


16QAM 


64QAM 


The proposed code 




16 


3.65 


6.20 


7.33 


The rate-5/4 punctured code in [12| 




12.8' 


4.81 


7.36 


8.49 



This value was obtained for 4-/16-QAM as stated in 1121 




-H — Proposed code, 4-QAM 

■ H - punctured P. Srinath e(.a/code, 4-QAM 
-V — Proposed code, 1 6-QAM 

■ V - punctured P. Srinath e(.a/code, 1 6-QAM 



6 8 10 12 14 16 16 

SNR per receive antenna 



Fig. 1. CER performance for 4x2 configuration and 4-/16-QAM 




-B — Proposed code, 4-QAM 

■ B - punctured P. Srinatti ef.a/code, 4-QAM 
-V — Proposed code, 1 6-QAM 

■ V - punctured P. Srinatti ef.a/code, 16-QAM 



6 8 10 12 14 
SNR per receive antenna 



Fig. 2. Average complexity for 4x2 configuration and 4-/1 6-QAM 

The coding gain S is equal to the minimum Coding Gain 
Distance (CGD) fp-J, or mathematically: 



S = min det 

s,s'ec 



(X(s)-X(s'))''(X(s)~X(s')) 



CGD(X(s),X(s')) 

inin |Jef[(X(As))]p 

AseAC\{0} 



(21) 



where As = s — s', AC is the vector space spanned by As. 
However, the code structure in ^TT) imposes: 

min Idet [(X(As))] I < Idet [0(2, 0, 0, 0, 0, 0] I = 16 

AseAC\{0} 

As a result, the angle <j) that maximizes the coding gain has 
to satisfy: 

|t/ef[(X(As))] I > 16,V As 7^ (22) 
For the proposed code we have: 

/ 6 \^ / 10 \ ^ 

\det[{X{As))]\J 

where Axi = 2ni 
/ 6 



\i=l 
10 



(23) 



, 10, and 



b = 2 



— 4|^(Ax7Aa;4 — AxgAx^ + AxqAxq ~ AxioAx2)^ 

+ (Ax7Aa;6 + Ax^Ax2 — Aa:9Aa;4 — AxiqAxs)^ 

+ (Aa;7Aa;2 — AxgAxg — AxgAxa + Aa;ioAx4)^ . 
A simplification of the expression of h is possible by noting 



that: 




Axl + Axl 



AxT^ 



Ax^ 



8 



Axi 



Ax^ 



[Ax^ 



L,g^ (Aa;^ - 

^1 -r AxX) {Ax^i + Aa;|) + {Ax\ + Ax^) {Axi + ^^^lo) 
Applying the Euler's four square on the first term and the 
Fibonacci's two square identities on the rest of terms we 
obtain: 

' 10 




VJ=7 

AxgAxa + Aa;9Aa;6 — Aa;ioAa;2)^ 
Aa;8Ax2 — Ax9Aa;4 — AxioAxa)^ 
+ (Aa;7Aa;2 — Aa;8Aa;g — AxgAxa + Aa;ioAa::4)^ 
+ (Aa;7Aa;3 + Aa;gAx4 + AxgAa;2 + AxioAxg)^ 
+ (Aa;7Aa;i + Ax^Ax^f + (A.tsAxi - Aa;7Aa;5)^ 
-h (AxgAxi + AxioAxs)^ + (AxioAxi - AxgAxs)^ 



i=l 



Therefore, one may write b in more compact form as below: 



2(a? 




(24) 



Setting X = e-'^'^, the discriminant A of the second degree 
equation ( |23] l is expressed as: 



A 





Consequently, the roots of equation ( [23| ) are: 

-b±j^ 



Ai, 



62 



Ax: 



A2 = At, |Ai| = |A2| = 



Ell A^f 



For the sake of simplicity, we will denote hereafter 
Y^^=i A2^| = CTi and X]j=7 Aa;^ — (T2. In the case of di 7^ 0-2, 
equation. [23] can be lower bounded as below: 



\det[(X{As))] 



As^O 

CTi ^(T2 



{x - Ai) (a; - A2) 

(N|-|Ai|)(|x|- 
2 



IA2I) 



(27) 



1 - 



^2 
\2 



= ((72 - (Tl) > 16 

where the latter inequality follows by substituting Axi — 
2ni, Hi ^ Z as we are dealing with unnormalized QAM 
constellations. If cri = (T2 = cr, equation ( (23] l can be written 
as 



As^O 



\2a^ cos(20) + 6| 



(28) 



al)). Taking cos(20) = 1/5, 



\det mAs))]\ 

where 6 = 2 (cr^ _ 2(af 
we have to prove that: 

2cr2 

+6 > 16 V CTi = (72 cr 7^ 0. (29) 

5 

Multiplying both sides by 5 and using ( |24| l, the above inequal- 
ity becomes: 



12(7^ - 20(a^ 



•i' \{ai 



,08)^0 

I we may have: 



;;2 



Denoting = ^ and a 
1 125-2 - 20(a^ - 
The above inequality is satisfied iff: 



> 5 X 16 (30) 
>5. (31) 



,08)^0 



1352 



> 2. 



(32) 



-2 -r-3;|(a,,...,a8)/0 

However, the L.H.S of the above inequality can be considered 
as a special case of: 



(33) 



(Xl,X2,X3,X4)#0 

where Xi £ Z. This type of equations have been extensively 
studied in the mathematical literature dealing with the solv- 
ability of quadratic Diophantine equations (see |15|). Applying 
theorem 6 in ifTSl we have: 

3Xl-5{Xi+X: 



'Xi)^Oy {Xi,X2,X3,Xi)^0 (34) 



1 (mod 8) with 3-5-5-5 = 4 (mod ; 



as —3x5x5x5 
Moreover, 

3Xf - 5(X| + xi + Xl) ^ ±1. (35) 
Otherwise, we must have: 

iXf EE ±1 (mod 5) (36) 
which cannot be true, since the quadratic residues modulo 
5 ai-e 0,1 and 4 [161, thus iXj = 0, ±3 or ± 2 (mod 5). 
Therefore, we can write: 



ixl - b{xl + xl + Xl) 

which in turns implies: 

|c/£'f[(X(As))] 

(26) thus completing the proof. 



(Xi,X2,X3,X4)#0 

> 16 

As^O 



> 2 



(37) 
(38) 
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